arXiv:l 502.07914v5 [quant-ph] 4 Aug 2015 


Invertibility of retarded response functions for Laplace transformable potentials: 
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A theorem for the invertibility of arbitrary response functions is presented under the following 
conditions: the time-dependence of the potentials should be Laplace transformable and the initial 
state should be a ground state, though it might be degenerate. This theorem provides a rigorous 
foundation for all density-functional-like theories in the time-dependent linear response regime. 
Especially for time-dependent one-body reduced density matrix (1RDM) functional theory this is 
an important step forward, since a solid foundation has currently been lacking. The theorem is 
equally valid for static response functions in the non-degenerate case, so can be used to characterize 
the uniqueness of the potential in the ground state version of the corresponding density-functional¬ 
like theory. Such a classification of the uniqueness of the non-local potential in ground state 1RDM 
functional theory has been lacking for decades. With the aid of presented invertibility theorem 
presented here, a complete classification of the non-uniqueness of the non-local potential in 1RDM 
functional theory can be given for the first time. 


I. INTRODUCTION 

The main statement of the Runge-Gross theorem is 
that the potential to density mapping is invertible mod¬ 
ulo a time-dependent constant in the potential jl]. This 
statement implies that knowledge of the time-dependent 
density and the initial state is in principle sufficient 
to fully characterize the evolution of a quantum sys¬ 
tem. The invertibility theorem by Runge-Gross there¬ 
fore forms the cornerstone of time-dependent density 
functional theory (TDDFT) (2}|5|- Unfortunately, the 
Runge-Gross theorem only holds for potentials which are 
Taylor-expandable in time. Taylor expandability of the 
potential is actually a too stringent condition and can be 
loosened as demonstrated by the invertibility theorem for 
linear response by Van Leeuwen (B| and more recently, by 
work of Tokatly on lattice systems IMj and the fixed- 
point approach by Ruggenthaler and Van Leeuwen na¬ 
ns]. 

Practical TDDFT calculations are almost exclusively 
performed with the help of an auxiliary non-interacting 
reference system, the Kohn-Sham system due which 
has the same density as the fully interacting system 
thanks to the exchange-correlation potential. The ex¬ 
act exchange-correlation potential depends on all densi¬ 
ties at earlier times and on the initial states in a com¬ 
plicated manner [14]. This history dependence of the 
exact exchange-correlation potential is neglected in prac¬ 
tice and replaced by its ground state version in which 
the instantaneous density is inserted. This neglect of 
memory dependence is known as the adiabatic approxi¬ 
mation and the results are typically very satisfactory for 
polarizabilities and local valence excitations. Especially 
when a good model for the exchange-correlation poten¬ 
tial is used even Rydberg excitations can be reproduced 
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reliably [T5HT7] . Practical TDDFT calculations fail for 
more complicated excitations such as charge transfer ex¬ 
citations UHl HI] and bound excitons miE], when the 
hole and electron are not localized close to each other [22], 
though some progress has been reported in their TDDFT 
description 12.'11125 1 . Even more problematic are dou¬ 
ble [26l [27] and bond-breaking (28[ [29] excitations. The 
main problem is that the density is not a natural quan¬ 
tity to describe these excitation processes and the non¬ 
interacting Kohn-Sham system is also of no avail. 

A more natural quantity to deal with these more com¬ 
plicated physical processes is the one-body reduced den¬ 
sity matrix (1RDM). In particular its fractional occupa¬ 
tion numbers are good descriptors of correlation effects. 
Indeed, it has been demonstrated that time-dependent 
1RDM functional theory is capable of correctly describ¬ 
ing charge-transfer excitations, double excitations and 
bond-breaking excitations [5PH52] even within the adi¬ 
abatic approximation. Unfortunately, no proper formal 
justification for time-dependent 1RDM functional theory 
has yet been published. The main purpose of this paper is 
to partially eliminate this caveat by presenting an invert¬ 
ibility theorem for non-local potentials and the 1RDM in 
the linear response regime. 

In previous work |3d| one tried to avoid the lack of 
a formal foundation by invoking the the Runge-Gross 
theorem. If all observables are functionals of the time- 
dependent density, they certainly are of the 1RDM, since 
the density can readily be extracted from the 1RDM. The 
problem is that there are (infinitely) many lRDMs gen¬ 
erating the same density, but only one of these lRDMs 
corresponds to the local potential belonging to that den¬ 
sity, the ‘local’ 1RDM. So if we want to use the Runge- 
Gross theorem (or one of the extensions) as a foundation 
for time-dependent 1RDM functional theory, we are only 
allowed to use these ‘local’ lRDMs. Since the character¬ 
ization of these ‘local’ lRDMs is seems to be impossible, 
direct use of the Runge-Gross theorem does not lead to 
any viable theory. 
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Actually, if we follow the philosophy of TDDFT more 
closely, we should not consider the mapping from local 
potentials to lRDMs, but from non-local potentials to 
lRDMs, since the non-local potential is the natural con¬ 
jugate variable of the 1RDM. With a non-local poten¬ 
tial I mean a one-body potential non-local in space, but 
still local in time. A generalization of the Runge-Gross 
theorem from local potentials and densities to non-local 
potentials and lRDMs would therefore be more appro¬ 
priate. Unfortunately, a straightforward generalization is 
not possible, since the commutator between the 1RDM 
and the interaction does not vanish, [ 7 , W\ 7 ^ 0. 

The invertibility theorem for the density response 
function for Laplace transformable potentials by Van 
Leeuwen [ 6 ] is much more amenable to generalization 
to the 1RDM. Assuming that the initial state is a non- 
generate ground state and that the time-dependent part 
of the potential is Laplace transformable, the theorem 
by Van Leeuwen states that the density response func¬ 
tion is invertible up to a constant shift in the potential. 
The proof by Van Leeuwen can be split into two parts. 
The first part consists of the derivation of a necessary 
and sufficient condition for a perturbation to yield zero 
response, i.e. necessary and sufficient condition for a po¬ 
tential to belong to the kernel of the response function. 
This first part does not need the special properties of the 
density operator and can therefore readily be general¬ 
ized to arbitrary operators, e.g. the 1RDM operator, the 
current operator, the (non-collinear) spin density opera¬ 
tor and the kinetic energy density operator. The second 
step of the proof by Van Leeuwen is to check the neces¬ 
sary and sufficient condition for the density operator, so 
the second part depends on the particular properties of 
the density. A generalization of the second step to arbi¬ 
trary operators is therefore not possible, but needs to be 
considered separately for each operator. In this article 
I will work out the second step for the 1RDM operator, 
which allows me to completely characterize the kernel of 
the 1RDM response function. 

The first part of the proof by Van Leeuwen requires at 
one point that the initial ground state is non-degenerate. 
This non-degeneracy requirement is quite a nuisance, 
since this requirement immediately excludes all open 
shell systems. I will not only repeat the first part of the 
proof by Van Leeuwen for arbitrary operators, but also 
show how the theorem can be extended to include degen¬ 
erate initial states. Including initial degenerate ground 
states leads to an additional condition which needs to 
be checked. I work out this additional condition for the 
density operator and show that the inclusion of degen¬ 
erate states does not lead to additional potentials in the 
kernel of the density response function. This provides 
an extension of the invertibility theorem for the density 
response function by Van Leeuwen to degenerate ground 
states. Degenerate ground states are also considered for 
the 1RDM response function. Including the possibility of 
degeneracy leads to a more general expression for some 
of the potentials in the kernel of the 1RDM response 


function found non-degenerate case. These results put 
linear response time-dependent 1RDM functional theory 
on a rigorous foundation. A more complete theoretical 
framework would be achieved by also addressing the v- 
representability question. The generalized invertibility 
theorem only determines to which extent the perturba¬ 
tions Sv yielding $7 are unique. It does not answer the 
question if there actually exists a Sv which yields the de¬ 
sired ($ 7 , i.e. the question of u-representability of £ 7 . For 
a complete theory, also the v-representability question 
should be addressed to specify exactly which 6 7 we are 
allowed to use in linear response time-dependent 1RDM 
functional theory. This paper focusses on the invertibility 
of the response function, so the u-representability ques¬ 
tion is beyond the scope of this article. 

As Van Leeuwen showed in Ref. [Mj, his invertibility 
proof for the density response function also works in the 
time-independent case by taking the limit of the variable 
in the Laplace transform to zero. Especially for ground 
state 1RDM functional theory this is an important re¬ 
sult, since a Hohenberg-Kohn like proof for the non-local 
potential to 1RDM mapping does not exist as pointed 
out by Gilbert [35]. Using the invertibility proof for the 
1RDM response function for the time-independent case I 
can give a full classification of the non-uniqueness of the 
non-local potential featured in 1RDM functional theory 
for the first time. Application to the time-independent 
case only works for non-degenerate ground states unfor¬ 
tunately, since degeneracies are treated in a fundamen¬ 
tally different manner in both cases. 

The paper is organized as follows. I start by repeating 
the first part of the invertibility proof by Van Leeuwen 
for arbitrary operators and extend his approach to han¬ 
dle degenerate ground states as well. After necessary and 
sufficient conditions have been obtained to characterize 
the potentials which do not lead to a response, I work 
these conditions out for the density and the 1RDM op¬ 
erators. First I consider the density operator to check 
against the result by Van Leeuwen [33J and to extend 
his result for the density response function to degener¬ 
ate ground states. Next, I consider the 1RDM operator 
to find the non-local potentials which do not lead to a 
response of the 1RDM. First I consider the simpler case 
of non-degenerate ground states and then take the ad¬ 
ditional necessary condition into account to extend the 
result to degenerate ground states. In the last part before 
concluding, I discuss the implications of these results for 
ground state 1RDM functional theory and point out in 
more detail why degeneracies need a different treatment 
in the ground state case. 


II. THE GENERALIZED INVERTIBILITY 
THEOREM 

Now let us repeat the first part of the invertibility proof 
by Van Leeuwen [5] for arbitrary operators Qi and gen¬ 
eralize the proof to degenerate ground states. The set of 
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operators {Q,} can be any set of self-adjoint operators of 
interest, e.g. dipole and quadrupole operators. The index 
i is also allowed to be a continuous index to represent a 
self-adjoint operator density such as the density opera¬ 
tor, n(r). Of course, a mixture of continuous and discrete 
indices is also allowed such as in the spin-density opera¬ 
tor, n(x), or the 1RDM operator y(x, x'), where x := r<r 
is a combined space and spin coordinate. I will use the 
symbol := throughout the paper to emphasize definitions. 

Having selected some set of self-adjoint operators and/ 
or operator densities of interest, we consider perturba¬ 
tions by these operator Qj with strengths 5vj(t'). Note 
that we need Svj(t') £ R to ensure that the total per¬ 
turbation remains hermitian. Now we consider the linear 
response of the expectation values of the same set of op¬ 
erators [SB] 

= Y ~] [ d ~ t')6vj(t'), ( 1 ) 

3 J ° 

where Xij (t ~ t') is the retarded/causal linear response 
function which can be defined as 

Xij(t - t') := -i 6{t - t')('bo\[QH 0 ,i(t),Q H o,j(t')]\'bo). 

( 2 ) 

In its definition we have used the operators in their 
Heisenberg representation with respect to the unper¬ 
turbed Hamiltonian, QH 0 ,i(t ) := e lHot Qie~ lHot , which is 
often referred to as the interaction picture. We have also 
used the Heaviside function, which is defined as 


6{x) 


1 for x > 0 
0 for x < 0 . 


The retarded response function can alternatively be ex¬ 
pressed as a sum-over-states (its Lehmann representa¬ 
tion) as [34ll36ll37] 

=ie(t-t')^2e inKit ~ t ' ) qf*qf : + c.c., (3) 

K 

where fix '■= Ex — E$ > 0 are excitation energies and 
fix = 0 only for K < D, so D denotes the multiplicity of 
the ground state degeneracy [3S]. Further, we have de¬ 
fined qf := (tE'olQil'I'A')- Note that the initial state can 
be excluded from the sum, since = (^olQil^o) G 
because the operators Qi should be hermitian. Inserting 
the sum-over-state expression for the response function 
in 01 the response of the expectation value of the oper¬ 
ator Qi(t) can now be written as 


6Qi(t) = i£ 9 f / d t'a K (t')e inK ^ () -t-c.c., 

K Jo 

where we have defined 


The integral has the form of a convolution product over 
the interval [0, t\ and can be transformed into a normal 
product by taking the Laplace transform 

E[6Qi](s) = i £ q** £ [ a ^o 3 > + c-c-, 

K S — 1“A 

where the Laplace transform is defined as 

pOO 

C[f](s) := / dte~ st f(t). 

Jo 

Now we multiply this equation by the Laplace transform 
of the potential £[(5ui](s) and sum over the index i to 
obtain 

= - 2 £ -^%^\C[a K \( S )\ 2 . 

t K S + S Lr 

In absence of response, we have that SQi = 0, so we also 
have that C[5Qi] = 0 and we obtain from the previous 
equation that for zero response we necessarily have 

° = E^TSr |£[wlW|2 ’ 

K b -T 1 l k 

Because Qk > 0 and only for K < D we have fix = 0, all 
the contributions for K > D are positive. Therefore, one 
necessarily has £[ok-](s) = 0 for K > D, so ax{t) = 0 
for K > D as well. Strictly speaking, we should say that 
ax(t) = 0 almost everywhere for K > D, since ax(t) ^ 0 
on a set of measure zero in time would not contribute to 
the integral of the Laplace transform. From its defini¬ 
tion Q it is clear that there are only two possibilities for 
ax{t) = 0 almost everywhere. The first possibility is the 
absence of a perturbation, 5vj(t) = 0 almost everywhere. 
If we further assume that we are only interested in the 
classical solutions of the Schrodinger equation (these are 
the usual physical wave functions defined at each point in 
time), we have as an additional condition that the poten¬ 
tial needs to be continuous up to its first order derivative 
in time, dvj (t) £ C 1 [SH], so ‘almost everywhere’ could be 
dropped. For more details on the solvability of the time- 
dependent Schrodinger equation, I refer the reader to an 
excellent introduction in Ref. m ■ This first possibility is 
trivial, and will be excluded from further discussion. The 
other possibility is that there exists one or more linear 
combinations of the operators under consideration 

L n = 

3 

such that ('F 0 |i/ n |'F A -) = 0 for all K > D. Note that the 
linear combination should remain hermitian, so 5v™ £ R. 
This implies that such a linear combination acting on the 
initial state, L n | H/q) , should not produce any components 
outside the degenerate subspace, i.e. 


^Kit) :=^qf5vj{t). 
j 


Ln 1 * 0 }= £ K), 

K<D 


(4) 


(5) 
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where 1% '■= ('F 0 \L n \ ^k) . In the case of a non-degenerate 
ground state the situation simplifies to an eigenvalue con¬ 
dition 

£„|*o> = i n |*o>- (6) 

In words, for a non-degenerate initial ground state, the 
response can only be zero nontrivially, if there exists a 
linear combination of the operators Qj for which the ini¬ 
tial state is an eigenstate. Note |\Ho) being an eigenstate 
of L n is sufficient, though not necessary for degenerate 
ground states, since L n |^o) is still allowed to have com¬ 
ponents in the degenerate subspace ([5]). 

Though we have shown that aKit) = 0 for K > D is 
necessary for absence of response, we also need to check 
if this condition is sufficient. Now suppose that condi¬ 
tion ([5| holds for some initial state |^o) and some oper¬ 
ator Q n . Note that we can always make a linear trans¬ 
formation of the set of operators such that the operators 
L n are all explicitly contained in the set {Q n }- In that 
case only the states with Q K = 0 will contribute to the 
sum-over-state expression in ([3]), so reduces to 

Xni{t-t') =i 0(t-t')^23m[q^*q^] +c.c. 

K 

= id{t - t'){^o\[Qi,Qn]\^o}- 

so as an additional requirement for zero response apart 
from = 0 for I\ > D, we find that 

<M&,&.]|*0> = 0 v,. (7) 

A number of remarks on this condition are in order. If 
the initial state is an eigenstate of the operator Q n , this 
condition is automatically satisfied. Hence, it is sufficient 
for non-degenerate ground states to check condition (| 6 |. 
This implies that only in the case of a degenerate initial 
state for which Qnl’I'o) has some components in the de¬ 
generate subspace, condition 0 needs to considered ex¬ 
plicitly. Since operators commute with themselves, con¬ 
dition ([T]) is trivially satisfied for i = n, so the check only 
needs to be performed for i 7 ^ n. 

Example. To get a feeling how condition ([T]) comes 
into play, consider a three dimensional Hilbert space, 
H = {|2s), \2p x ), \2p y )}. We take the Hamiltonian of 
the hydrogen atom as our initial Hamiltonian, Hq, so the 
states in H are degenerate. The condition 0 is therefore 
trivially satisfied and only the additional condition Q 
needs to be considered. For the operators Qj : H — > H 
we consider the unit operator and the dipole operators 
in the x- and y-direction, {Qj} = {l,x, y}. We select 
the 2s-orbital to be our initial state, |^o) = |2s). Since 
any state is an eigenstate of the unit operator, we im¬ 
mediately find that at least the unit operator is a part 
of the kernel of the response function. The dipole opera¬ 
tors produce an additional component in the degenerate 


subspace gO] 

x\2s) ■= = —3|2p x ), 

tlnen 

y| 2 s) := ^2\ijji){'ip i \y\2s) = —3|2 p y ), 
linen 

so the additional condition ([TJ) needs to be checked ex¬ 
plicitly. The dipole operators already commute among 
themselves, [x,y\ = 0 , so also the dipole operators be¬ 
long to the kernel of the response function. Because all 
operators belong to the kernel of the response function, 
we actually have \ = 0- Such a response function is 
rarely encountered in practice and is a consequence of 
the special choice of the Hilbert space, H, and the set of 
operators {l,x,y}. If, for example, also the |3 p x ) state 
would be included in the Hilbert space, the x-operator 
would produce also components outside the degenerate 
subspace 

97 (S48 

x\2s) = ^^\ipi){ilJi\x\2s) = ^ 5625 l 3 fe) - 312^). 

ilnenu{\3 Pn: )} 

By the first condition ([ 5 ]), the x-operator would not be¬ 
long to the kernel of the response function anymore. 
Note that the extension has only consequences for the 
x-operator, so the unit-operator and y -operator remain 
in the kernel of the response function. 

The x-operator can also be lifted out of the kernel 
of the response function by adding additional operators. 
For example, consider an extension of the set of opera¬ 
tors by an operator which applies the momentum oper¬ 
ator three times in the x-direction, id}}. = (— \d x ) 3 m- 
The action of id% on the initial state is defined to be 

id 3 |2s) := ^|Vk)(^|i<9 3 |2s) = ^-\2p x ). 
ipie n 

Consider now the additional condition 0 for id Only 
the dipole operator in the x-direction yields a non¬ 
vanishing commutator and its expectation value for the 
initial state gives 

<2s|[z,i3 3 ] l 2s ) = °- 

The operator id^. is therefore not part of the kernel of the 
‘extended’ response function. Note that this result also 
implies that for the dipole operator in the x-direction 
condition Q is not satisfied anymore. Only the unit 
operator and the dipole operator in the y-direction span 
therefore the kernel of the ‘extended’ response function. 

III. DENSITY RESPONSE 

As a minor check, let us consider the density response 
function to see if we recover the original result by Van 
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Leeuwen [B], For the density response function our op¬ 
erators are Q r = n(r), where n(r) := '0' (ra)i/j(ra). 

The only linear combination for which a non-degenerate 
ground state is an eigenstate is the number operator 

N := Jdrh(r). (8) 

Only if the density would vanish in some region for |To), 
there would be other linear combinations for which |dr 0 ) 
would be an eigenstate. This possibility is typically ex¬ 
cluded in DFT so we recover the same result as 

Van Leeuwen [B] that only a spatially constant potential 
gives a zero density response. 

Now let us investigate the consequences of a degener¬ 
ate ground state by considering only one particle first. A 
non-constant potential yielding a zero density response 
is readily constructed as vk{ r) = Ta (r)/To(r) for 0 < 
K < D, which by construction satisfies <©• This con¬ 
struction only works if the wave function d'A'(r) has no 
imaginary part. Otherwise the potential vk would not 
be hermitian. Assuming that the ground states 'FaM 
are indeed real, the initial state should at least have one 
nodal surface to allow for the degeneracy. Further, be¬ 
cause the states need to be orthogonal, not all of their 
nodal surfaces should coincide. The potential Vk would 
therefore be infinite along some nodal surface of T 0 {TBj . 

For the next step we need to take additional conditions 
on the allowed perturbation into account, to ensure that 
the Hamiltonian remains self-adjoint p2i !47j . For po¬ 
tentials over M 3 one requires the perturbations to be in 
the class of Kato potentials K. := L 2 (R 3 ) + L°°(K 3 ) (15) . 
This means that the perturbing potential needs to be de¬ 
composable in parts which are either bound or square 
(Lebesgue) integrable. For example the Coulomb poten¬ 
tial can be split into two parts as 

1 0(1-|r|) fl(|r|-l) 

M M |r| 

The first part contains the Coulomb singularity, which 
that is square-integrable, so the first part belongs to L 2 . 
The second part is the outer region of the Coulomb po¬ 
tential which is bounded (between 0 and 1), so belongs 
to L°°. The Coulomb potential is therefore a proper per¬ 
turbative potential, since it is a Kato perturbation. An 
example of a potential which is not Kato is the har¬ 
monic oscillator potential, ^uj 2 \r\ 2 . Due to its diver¬ 
gence for | r| —> oo the outer part of the harmonic po¬ 
tential is neither bound nor square integrable. Return¬ 
ing to our potentials vk{ r) = \E , A'(r)/^ f o( r ), we expect 
these potentials to behave as 1 /z in the direction orthog¬ 
onal to the nodal surface. These potentials are there¬ 
fore expected not to be square integrable in regions in¬ 
cluding some part of the nodal surface. To proof this 
suspicion, we use the result by Kato [48] that the solu¬ 
tions of the time-independent Schrodinger equation are 
continuous over whole K 3JV and their derivatives locally 
in L°°, i.e. the solutions are locally Lipschitz. Locally 


Lipschitz means that for each xo,yo £ R 3JV there ex¬ 
ists some neighborhood V and a constant Ky such that 
l’F(x) — T(y)| < AV|x —y| for any x,y £ V. Now take 
a neighborhood around some point at the nodal surface 
of 4>o(r) and call this neighborhood S. The local Lips¬ 
chitz condition for To simplifies in this neighborhood to 
|To(s)| < Ks\s\ where s denotes the distance from the 
nodal surface. The sign of s indicates on which side of 
the surface we are. Since the nodal surfaces of T/f(r) and 
To(r) do not coincide, we can always choose some neigh¬ 
borhood such that min(|TA-|) = E > 0. The L 2 -norm of 
the potential within S can now be estimated as 


dr 


^g(r) 

4'o(r) 


> 




^ area(5'*) 



where area(5*) > 0 and e > 0. In the last inequality, S* 
is a part of the nodal surface contained in S, such that 
S* x [—e, e] C S. Thus the last integral is basically over 
a slab of thickness 2e contained in S along the nodal sur¬ 
face. The inequality therefore shows that the potential 
vk { r ) is not square integrable in the region S, so this 
part of the potential is not in L 2 . Since the potential 
is obviously not bounded in this region, i.e. not in L °°, 
the potential vk{ r) is not a Kato perturbations. The 
potentials vk{ r) are therefore not admissible, so degen¬ 
erate initial states do not form any complication for the 
invertibility of the density response function. The final 
result is that also in the degenerate case the kernel of the 
density response function only consists of a constant po¬ 
tential. It is obvious that the same conclusion also holds 
for more than one particle. 


IV. 1RDM RESPONSE 

We will now consider the 1RDM response function. 
The 1RDM operator is defined as 

7 (x,x / ) := ^ t (x , )V>(x), 

where i/;(x) and ^(x) are the usual field operators and 
x := r a is a combined space-spin coordinate. We will first 
limit ourselves to a non-degenerate ground state as initial 
state, since this case already leads to several situations 
which need to be considered. 


A. Non-degenerate ground state as initial state 

Because the density is simply the diagonal of the 
1RDM, n(r) = 7( r<7 i rcr )> the constant potential is 

also present in the kernel of the 1RDM response func¬ 
tion. Since the 1RDM contains more flexibility than the 
density, one would expect that there are more possible 
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potentials that give a zero response than only the spa¬ 
tially constant potential. Indeed, any one-body operator 
can be represented by the 1RDM, so if the initial state is 
an eigenfunction of some one-body operator, this oper¬ 
ator is also present in the kernel of the 1RDM response 
function. 

In particular, the non-relativistic Hamiltonian does not 
depend on spin, so a non-degenerate ground state is nec¬ 
essarily a singlet state. This implies that the ground 
state is an eigenstate of the total spin-projection opera¬ 
tor in arbitrary directions, Sj\&o) = 0- Since the total 
spin-projection operator can be expressed as a one-body 
operator, it is also part of the kernel of the 1RDM re¬ 
sponse function. Note that this situation also occurs in 
spin-DFT [15115U] . 

Since symmetry in the system implies that the Hamil¬ 
tonian commutes with one or more symmetry opera¬ 
tors, the eigenstates of the Hamiltonian can be chosen 
to be eigenstates of some of those symmetry operators as 
well. Therefore, one would expect that also these sym¬ 
metry operators belong to the kernel of the 1RDM re¬ 
sponse function. However, the Coulomb interaction of 
the Hamiltonian couples all the particles, so these sym¬ 
metry operators need to be many-body operators in gen¬ 
eral. Continuous symmetries form an exception, since 
their generators can be expressed as one-body operators. 
For linear molecules this would be the rotation around 
the z-axis, i.e. the L z operator. Atoms would also in¬ 
clude the other total angular momentum operators, L x 
and L y . For systems which are homogeneous in one or 
more directions, e.g. the homogeneous electron gas, the 
corresponding momentum operator(s) would also be part 
of the kernel of the 1RDM response function. 

To proceed with the analysis, we will work in the natu¬ 
ral orbital (NO) basis of the 1RDM of the initial ground 
state, which can be obtained by diagonalizing the 1RDM 

7(x, x') = n k </>fc(x)</>i((x'). 

k 

The eigenvalues are called the (natural) occupation num¬ 
bers and the eigenfunctions are the natural orbitals m- 
The occupation numbers sum to the total number of 
particles in the system, TV, and for fermions they obey 
0 < n fc < 1. The integer values are special, since 
rik = 0 implies that the NO </>fc(x) is not present in 
any determinant in the expansion of the wavefunction, 
aJ*o) = 0 , where a k is the annihilation operator for the 
NO </>fc(x). Likewise, a fully occupied NO, nk = 1, im¬ 
plies that the NO ^fe(x) is present in all determinants, so 
4 a fc ]*o) = |*o> M- From these properties, we readily 
find that 


7iul*o) 


'01*} 

< l|tt> 

Y Ck ^ k ) 

■ K 


if nk = 0 or 
ri; = 1 and k ^ l 
if 77 ; = 1 and k = l 
otherwise, 


where the 1RDM operator is now represented in the NO 
basis, jk,i : = a k . Hence we find that the ground state is 
an eigenstate of the 1RDM operator if ti^ = 0 or ri; = 1. 
However, we have to keep in mind that the potential 
should be hermitian, so if 5vu / 0, also 8v* k ^ 0. Thus 
for the state |*) to be an eigenstate of both 7 kj i and 
7z,fc, we additionally need that nk = 1 or m = 0. This 
situation can only occur if nk = ni = 0 or nk = n; = 
1. We find therefore, that the perturbations within the 
fully occupied block or within the completely unoccupied 
block have a zero response in the 1RDM, as is actually 
well known for non-interacting systems, e.g. the Kohn- 
Sham system in DFT. Note that this discussion includes 
the one-particle case, since that is also non-interacting. 

For interacting systems, the occupation numbers are 
predominantly fractional, 0 < nk < 1, and for Coulomb 
systems there is strong evidence that they all are 17?- 
155] . One would expect that another special situation can 
occur if these fractional occupation numbers are degen¬ 
erate. To investigate this situation, consider the NOs as 
a basis and assume that </>i(x) and (02 (x) are two degen¬ 
erate NOs. The contribution of these degenerate NOs 
to the initial state can be made explicit by writing the 
initial state as 

l*o) = a\al\$%_ 2 ) + 4l*]v-i> + 4+ |*at>, 

where = 0 for i = 1 , 2 and any b £ {0, 1 , 2 , 12 }. 

Note that the states y 8> a M ) are not normalized in general. 
The action of the lRDM-operator on the initial state can 
be worked out as 


7ij|*o) = aj|*Ar-i) + 

for i,j = 1,2. Since the |*at) component vanishes, the 
only way that |'F 0 ) can be an eigenstate is to have the 
eigenvalue zero, so all components |*^) need to be can¬ 
celled. The components aj|*) v _ 1 ) are not present in the 
initial state |*o)- This follows from the fact that 
and (£>2 (x) are NOs, so 712 = 0. Since one generally can 
not eliminate these components by taking linear combi¬ 
nations of 7 ij, fractional occupation degeneracies do not 
cause additional potentials in the kernel of the 1RDM 
response function in general. 

A special situation occurs if |*) v _ 1 ) = 0. The only 
known interacting case is the two-electron system. The 
two-electron state in the NO representation can be writ¬ 
ten as an expansion of NO pairs to which each NO con¬ 
tributes only once [561459] 


|*o) =^2cka[al\). 

k= 1 

The coefficients in the expansion are called natural am¬ 
plitudes and are related to the occupation numbers as 
|c fc | 2 = Tik = n k . In the case of a singlet state, the NO 
pairs only differ in their spin part, (f>k( x ) = cf>k(r)a(a) 
and (f>k( x ) = <( , fe(r)/T(o’). In the triplet case the NO pairs 
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have different spatial parts and their spin parts can be 
identical. The paired NOs are degenerate and since we 
now have = 0, we find that 

° = 7fc,fcl^o) =7fc, fc |^o) = (%,k -7fc,fe)l^o), (9) 

so perturbations with these operators yield a zero 1RDM 
response [60 . 

The special structure of the two-electron state also 
causes other NOs with degenerate natural occupation 
numbers to yield zero response. For example, consider 
the contribution of two pairs of NOs to the initial state 

CiaJdjI) + C2Ci\d\\). 

Now we work out the action of the following perturba¬ 
tions which mix these NO pairs 


B. Including degeneracies 

Now let us consider if additional potentials will be part 
of the kernel of the 1RDM response function if we allow 
for degenerate ground states. We will do this by first 
assuming that there exists some non-local one-body po¬ 
tential which only creates components in the degenerate 
subspace when acting on the initial state 

#l*o> = ^ °’ ( 12 ) 

i,j 0 <K<D 

where uk = (’F 0 |f7|’I , if). Subsequently we check whether 
the additional necessary condition 0 is satisfied. This 
condition attains the following simple form in the NO 
basis 


(« 2 i 7 i ,2 +^2172,1)1^0) = V 2 \ c \ a \ a \\) + |), 

(* 7272,1 +^I 27 i, 2 )I^o} = *72 C 2 <i\<i\\) + vf 2 Cidja^l). 

Note that we added a second term to each perturba¬ 
tion, to ensure that the operator is hermitian. We see 
that both perturbations produce exactly the same deter¬ 
minants, so by combining both perturbations we might 
be able to cancel both. To eliminate the first deter¬ 
minant, aJjdjI), we need to set v \ 2 = — V21C1/C2 . To 
eliminate the second determinant, d|dg|), we need to set 
u^2 = — *’2i c 2/ci- This only works when the natural oc¬ 
cupations are degenerate, n\ = |ci| 2 = |c 2 1 2 = 112- In 
that case the following potential belongs to the kernel of 
the 1 RDM response function 

*>2i (71,2 - e iai 2 72 ,i) + v* 21 (72,1 - e _iai2 7i i 2), (10) 

which depends on the relative phase of the natural ampli¬ 
tudes, e‘“ 12 := C1/C2. It is readily checked that degener¬ 
acy between the NO pairs implies that also the potential 

*'21 (712 + e lai 2 72 l) + v 2i (721 + e lai 2 7 i 2 ) ( 11 ) 


belongs to the kernel of the response function. Note that 
the relative phase of the natural amplitude is important 
in the construction of the potentials (10) and Since 
the natural amplitudes can only be defined for a two- 
electron system, the notion of relative phases only makes 
sense for two-electron systems. The relative phases there¬ 
fore emphasize the special status of the interacting two- 
electron system concerning NO degeneracies. 

All non-local one-body potentials in the kernel of the 
1RDM response function have now been characterized for 
both the non-interacting case and the fully interacting 
Coulomb system. However, if one works in a small finite 
basis or uses some effective interaction which only af¬ 
fects some subspace, some special structure in the ground 
state might arise which causes additional potentials to be 
present in the kernel of the 1RDM response function. A 
complete proof including these cases would therefore re¬ 
quire additional assumptions or a more extensive analysis 
which would depend on the specific details. 


0 = ('Fo| [iki, U] I’Fo) = (rij - n k )u k i Vk,i- (13) 

This is a very interesting expression, since it tells us that 
only potentials U which have only non-zero matrix ele¬ 
ments coupling degenerate NOs yield a zero 1RDM re¬ 
sponse. This is a very stringent condition, especially in 
combination with the requirement that t/|H/o) is only al¬ 
lowed to have components in the degenerate subspace of 
ground states. 

Now let us investigate a number of systems of inter¬ 
est. First consider a system of non-interacting parti¬ 
cles with degenerate ground states. These ground states 
can be constructed by first solving effective one-particle 
Schrodinger equations h\(j>i) = ei\4>i). Assuming that the 
orbital energies are ordered in increasing order, Ci < e»+i, 
the complete span of degenerate ground states can be 
constructed by assembling all determinant with the low¬ 
est orbital energies 

l^u,..,w- fc ) = a ir" a L- fc a l'" a il>’ 


where «i < i 2 < ■ ■ ■ < iN-k € V := {fc+l, k+ 2,..., k+d} 
and d denotes the number of degenerate orbitals with or¬ 
bital energy ek+i = £k +2 = ■ • • = tk+d.■ Note that this set 
of degenerate ground states is not unique. Arbitrary uni¬ 
tary transformations among the degenerate ground states 
yield different spans of the ground state subspace which 
are equally valid. In particular, the initial ground state 
of the response function can be chosen as any linear com¬ 
bination of the degenerate ground states. 

A potential U satisfying condition (12) is readily con¬ 
structed by setting Uj = 0 if any i, j £ V and choos¬ 
ing some u.ij 0 for both i,j £ V. for the poten¬ 
tial U to belong to the kernel of the 1RDM response 
function, we additionally need that n, = rij for all ele¬ 
ments Uij 7 ^ 0. Note that the cases rij = rij = 1 and 
rii = rij = 0 were already covered before in the non¬ 
degenerate case, since |\R 0 ) will be actually an eigenstate 
of the potential U. The treatment of degenerate non¬ 
interacting ground states extends this result to any po¬ 
tential coupling only degenerate NOs. So all potentials 




with Uij = 0 for rii ^ rij will belong to the kernel of the 
1RDM response function in the non-interacting case. 

Now let us consider a system with a spin degenerate 
ground states, \S,M). These states are eigenfunctions 
of the S z operator, S Z \S, M) = M\S,M }, so the S z op¬ 
erator immediately belongs to kernel of the 1RDM re¬ 
sponse function. By operating with the S± operators 
we can obtain other states in the degenerate subspace, 
S±\S,M) = C±(S,M)\S,M ± 1). In second quantiza¬ 
tions, these raising and lowering operators can be ex¬ 
pressed as 

s+ = J2 “LA/? and = 4/A«- 

k k 


The operators S± are not hermitian operators, but we 
can make two independent hermitian combinations which 
are properly hermitian 


S X = ±(S+ + S-) 


and S v = i(S + -S_). 


Since these operators only produce components in the 
degenerate subspace, we have found proper potentials U 
as in (12). Now we need to check whether these opera¬ 


tors satisfy (13). We see that the S x and S y operators 


couple the different spin components of each spatial or¬ 
bital, so we need Uk a = nkp for all k for (131 to hold. 
This degeneracy only occurs for M = 0, so only in the 
case that S, 0) is the initial state, the S x and S y oper¬ 
ators also belong to the kernel of the 1RDM response 
function. Combined with our result for non-degenerate 
states, this means the S z is always part of the the ker¬ 
nel of the 1RDM response function for Hamiltonians not 
depending on spin. If additionally the system is spin- 
compensated, i.e. rika = n<k /3 for all k, the S x and S y 
operators are also part of the kernel, irrespective if the 
ground state is degenerate or not. Note that the same 
considerations also hold for the angular momentum oper¬ 
ators L if the Hamiltonian is invariant under all rotations, 
e.g. atoms and the homogeneous electron gas, though we 
need to check for different degeneracies in the occupation 
spectrum. For example consider an atom. The 3-axis can 
always be chosen such that the ground state is also an 
eigenstate of the L z operator. For the L x and L y op¬ 
erators to be part of the kernel of the 1RDM response 
function as well, we need nk t i >m = which implies 

that the 1RDM will be unperturbed when we make ro¬ 
tations around an arbitrary axis. 


C. Ground 1RDM functional theory 

The generalized invertibility theorem for the non¬ 
degenerate case is equally valid for the time-independent 
response function by taking the s —> 0 limit of the 
Laplace transformed quantities. The generalized invert¬ 
ibility theorem therefore provides the perfect opportunity 
to give a better classification of the uniqueness of the 


mapping from non-local one-body potentials to lRDMs, 
si-)). As Gilbert already mentioned in 1975 [35], the 
class of potentials which map to the same ground state 
1RDM will be larger than in DFT, but to the author’s 
knowledge no attempt has been made to give a full clas¬ 
sification of this non-uniqueness. We will show that the 
kernel of the time-dependent 1RDM response function 
exactly corresponds to the non-uniqueness of the non¬ 
local potential in ground 1RDM functional theory in the 
non-degenerate case. 

As Gilbert already showed [32], the second part of 
the Hohenberg-Kohn theorem can straightforwardly be 
generalized to lRDMs: the 1RDM of a non-degenerate 
ground state is unique. In other words, consider all the 
ground states corresponding to different non-local poten¬ 
tials, then there is a one-to-one correspondence between 
the non-degenerate ground states and their correspond¬ 
ing lRDMs. 

Now assume that there are two (non-local) poten¬ 
tials, V\ and # 2 , yielding the same non-degenerate ground 
state. Since the Schrodinger equation is linear, the po¬ 
tentials i>\ = (1 — A)fii + \v 2 yield exactly the same non¬ 
degenerate ground state. The set of potentials which 
yield the same non-degenerate ground state is therefore 
(simply) connected. To determine this set, it is there¬ 
fore sufficient to consider a perturbation to one of these 
potentials and check which potentials do not lead to a re¬ 
sponse to any order. As we have shown before, the first 
order 1RDM response only vanishes if the ground state is 
an eigenstate of the perturbation, but this also immedi¬ 
ately implies that the response will vanish to any order. 
We can therefore conclude that the kernel of the 1RDM 
response function exactly coincides with the class of po¬ 
tentials yielding the same ground state 1RDM. More pre¬ 
cisely, two non-local one-body potentials yield the same 
ground state (1RDM) if and only if their difference is 
part of the kernel of the 1RDM response function. Note 
that these considerations are not special for the 1RDM, 
but can be applied to any density-functional-like theory 
for which we are able to characterize the kernel of the 
response function. 

The degenerate case is beyond the scope of this ar¬ 
ticle. The main reason is that the degenerate case is 
handled in a fundamentally different manner in time- 
dependent and time-independent perturbation theory. 
Time-independent perturbation theory is based on the 
time-independent Schrodinger equation, which is an 
eigenvalue equation. The ground state is therefore only 
specified up to its degenerate subspace, from which an ap¬ 
propriate |T 0 ) needs to be chosen. This is reflected in its 
perturbation theory, since we need to diagonalize the per¬ 
turbation in the degenerate subspace and take the lowest 
eigenvalue. The perturbed state therefore depends on the 
direction of the perturbation as illustrated by the lowest 
energy surface in Fig. [T] Time-dependent perturbation 
theory is based on the time-dependent Schrodinger equa¬ 
tion, which is an initial value problem. The initial state 
|To) is therefore completely specified from the start, even 
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Figure 1. A sketch of the dependence of the lowest eigenstates 
of some quantum system as a function of a perturbation. The 
energy dependence of the initial state in the time-dependent 
case is shown by the thick line and is simply a linear func¬ 
tion. Time-independent perturbation theory always selects 
the lowest perturbed state, so I'&i) on the negative side and 
lit 2 ) on the positive site up to the next degeneracy point. 
To emphasize the jumps at the degeneracy points, the lowest 
energies are colored red. 

in the degenerate case. Therefore, contrary to the time- 
independent case, we do no need to (and can not) diag¬ 
onalize the perturbation in the degenerate subspace to 
select an appropriate zeroth order state, since it is sim¬ 
ply specified from the start as |d>o). The time-dependent 
situation is illustrated in Fig. [T] by the thick line. Due to 
the fundamental difference in dealing with degeneracies, 
the result for the time-dependent response function does 
not straightforwardly carry over to the time-independent 
response function in the degenerate case and a separate 
treatment is required. 

V. CONCLUSION 

To summarize, I have generalized the first step of the 
invertibility theorem for the density response function by 
Van Leeuwen [Bj, to arbitrary operators and to degener¬ 
ate ground states. For the nontrivial absence of response, 
it is sufficient that initial ground state is an eigenstate 
of the perturbation operator and also necessary in the 
case of a non-degenerate ground state. For a degenerate 
ground state, however, the action of the perturbation op¬ 
erator is allowed to yield additional components in the 
degenerate subspace, though the expectation value of the 
commutator of the perturbation with any operator under 
consideration needs to vanish to yield zero response as an 
additional condition 0. 

The theorem can be used to establish density- 


functional-like theories in the time-dependent linear re¬ 
sponse regime. The restriction to ground states is not 
very severe, since this is the initial state which is used al¬ 
most exclusively in practical linear response calculations. 
The determination of the kernel of the time-dependent re¬ 
sponse function also immediately carries over to the time- 
mdependent response function if the initial ground state 
is non-degenerate. This result is useful to establish time- 
independent density-functional-like theories. The kernel 
of the response function exactly coincides with the v t— > Q 
mapping in the non-degenerate case. The non-uniqueness 
of the potential can therefore be fully characterized as in 
the first Hohenberg-Kohn theorem for DFT. 

The generalized invertibility theorem has been applied 
to the density response function and it has been estab¬ 
lished that only the spatially constant potential belongs 
to its kernel, even for a degenerate ground state. Apply¬ 
ing the theorem to the 1RDM response function revealed 
that not only the constant time-dependent shift is part of 
the kernel, but also generators of continuous symmetries 
are possibly included. For non-relativistic Hamiltonians 
this would always be the S z operator and if the NOs are 
degenerate in both spin channels, rika = n/c/ 3 , also the 
other components of S belong to the kernel of the 1RDM 
response function, cf. non-collinear spin-DFT [ij5) f50] . 
Also the angular momentum operators are possibly in¬ 
cluded in the kernel of the 1RDM response function if the 
Hamiltonian is invariant under the corresponding rota¬ 
tions. The additional condition (J7| requires the relevant 
NOs to be degenerate as well. It is obvious that when 
spin-orbit coupling is included, the relevant operators to 
be considered would be J := L + S instead of L and S 
separately. For homogeneous systems, e.g. the homoge¬ 
neous electron gas, also the momentum operators —iV 
need to be considered. Further, the matrix elements of 
the non-local potential which couple within the fully un¬ 
occupied block or within the fully occupied block also do 
not lead to a first order response. In the non-interacting 
case actually all potentials coupling only degenerate NOs 
belong to the kernel of the 1RDM response function. Due 
to the intimate relation between a two-electron state and 
its 1RDM, degeneracies of the natural occupation num¬ 
bers give rise to additional non-local potentials in the 
kernel of the 1RDM response function, whose matrix el¬ 
ements couple the natural orbitals within the degenerate 
sub-block. This result not only puts time-dependent lin¬ 
ear response 1RDM functional theory on a rigorous basis 
but is also of high importance for ground state 1RDM 
functional theory, because it allows for a full character¬ 
ization of the non-uniqueness of the non-local potential 
for non-degenerate ground states for the first time. 
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